Abstract 



The spectral function of pions interacting with a gas of nucleons 
and A33-resonances is investigated using the formahsm of Thermo 
Field Dynamics. After a discussion of the zero A- width approximation 
at finite temperature, we take into account a constant width of the 
resonance. Apart from a full numerical calculation, we give analytical 
approximations to the pionic spectral function including such a width. 
They are found to be different from previous approximations, and 
require an increase of the effective A- width in hot compressed nuclear 
matter. The results are summarized in an effective dispersion relation 
for interacting pions. 
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1 Introduction 



One of the main experimental goals of modern nuclear physics is the investi- 
gation of hot, compressed nuclear matter (HCNM). Theoretical descriptions 
of HCNM are mostly based on, or at least motivated by, methods of rela- 
tivistic quantum field theory. However, within this framework it is a well 
established fact, that naive perturbation theory breaks down at finite tem- 
perature [jl], 0]. This is due the modification of space-time symmetry in the 
presence of matter or a heat bath, i.e. to the absence of stable asymptotic 
states for the observable physical particles To overcome this problem, 
first of all one has to use one of several existing formulations of quantum field 
theory with 2x2 matrix-valued propagators. 

While so far the Closed-Time-Path method (CTP) has been used to derive 
transport equations for heavy- ion collisions the model called Thermo 

Field Dynamics (TFD) offers technical advantages over CTP |^, 0]. Further- 
more, within such a formulation, proper asymptotic conditions necessary for 
a perturbative approach can be defined only in TFD (see ref. P for details). 

In this paper, TFD is applied to a system potentially interesting for the 
theoretical description of HCNM, namely to an interacting gas of nucleons, 
pions and Asa-resonances 0. This A-hole model is often used |]TU|, |TT|, |T^], 
but in our view the analytical structure of the pion propagator is far from 
well understood. 

We first describe the model as found in many applications, i.e. in the 
quasistatic zero-A-width approximation. Then a more realistic Lorentz-type 
spectral function for the A's is discussed, and the full pion propagator is 
calculated from a dispersion integral. 

To obtain a systematic extension of the simple A-hole model to finite 
A-width and finite temperature, analytical results are given for the pion 
propagator by introducing an asymptotic expansion of the A spectral func- 
tion to first order in its width. These results are simple enough to be used in 
further applications of the model, but still complete enough to contain the 
full information of the A-hole model. 

It is discussed, how to fit this simplified model to the full calculation, 
and how good the approximations are at various temperatures, densities and 
momenta. Finally, the model is used to extract an "effective" pion dispersion 
relation in HCNM. 

In the framework of TFD, the thermal instability of observable states can 
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be absorbed into a Bogoliubov transformation also for interacting systems. 
This Bogoliubov transformation defines stable quasi-particles which serve as 
basis for a perturbation expansion. It can be written in a single matrix form 
for the bosonic and the fermionic sector of the model as 

BsAn) = ( <'^i"' 7' ) . (1) 

The Bogoliubov parameters n are the phase-space distribution functions for 
bosons or fermions. While one has some freedom in parameterizing this 
transformation, the above form was found to be the most useful in ref. [0], 
since it makes the propagators linear in n. 

In the fermionic sector, i.e., for nucleons and A's, the n are Fermi-Dirac 
functions ^ 

with inverse temperature f3. The baryon number in each small volume and 
hence the baryon density is a constant parameter of the calculations. For- 
getting about the interaction at the moment, we express it in terms of bare 
"on-shell" energies 



En{p) = ^Jp^ + MI and E^{p) = + Ml , (3) 
and obtain as the baryon density without interaction 

= ^ l-^3MENip)) + lQ J-^n^iE^ip)) . (4) 

In the fermion propagators we neglect anti-particle states. The resulting 
2 X 2-matrices are 

4"a(Po,p) = -^ JdEAN,AiE,p) X 

/ I \ 



Po - E + ie 



1 



Bpin^^AiE)) (5) 



V po — E — ie / 

for both nucleons and A's. In this equation, An,a is the spectral function 
of nucleons or A's, and BF{nN,A) are Bogoliubov matrices containing the 
thermal information as defined above, = diag{l, —1). 



4 



2 Full pion propagator 

The full pion propagator is the solution of a Schwinger- Dyson equation, with 
the free propagator and a self energy function (polarization function) as in- 
put. Each of these quantities is 2 x 2-matrix valued, and the full as well as 
the free propagator obey 

It is easy to see, that this implies for the polarization function 

nil ^ ^ + = . (7) 

In the following, we use results derived in ref. P). From the linear relation (|^) 
between the matrix elements one obtains, that the matrix valued polarization 
function can be brought to triangular form by multiplication with bosonic 
Bogoliubov matrices Bb 



nii + ni2 (i + n,(E))ni2-n,(E)n2i 
nil ^ n^i 



(8) 



for any value of nT^{E). The diagonal elements are the retarded and advanced 
components of the polarization 

= n^^ + n^i . (9) 

As one can derive, the off-diagonal components of 11 are related by 

U^\E, k) - e-^(^-^^+'^^) U^\E, k) = Q. (10) 

Note, that this relation is quite similar (and indeed equivalent) to the Kubo- 
Martin-Schwinger boundary condition for propagators, which is equivalent 
to an equilibrium condition. 

The full pion propagator therefore can be diagonalized by Bogoliubov 
matrices, if and only if n.,^{E) is chosen such that the off-diagonal element in 
(H) vanishes, i.e. by setting 
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Since /i^r = /ia — /^at in our simple model, this is the pionic equilibrium 
distribution function, and thus the pion propagator is diagonal (separately 
for positive and negative energy states) only in case the pions have the same 
temperature as the baryons. For the results presented in this paper we have 
used = 0, but the equations can be applied also to cases with nonzero 
pion chemical potential. 

The full pion propagator is 



(A;o,fc) = JdEA^{E,k) X 



1 



V 



kn — E + ie 



\ 



E 



It 



^ BB{n^{E)){t^2) 

E-it J 

\ \ 



V 



kQ + E + ie I 



J 



A-„{E, k) is a positive function, the spectral function of the pion field. Thus, 
apart from the task of determining this function, (12) already solves the 
problem of pion propagation in HCNM. 

The spectral function is related to the imaginary part of the retarded and 
advanced propagator, 

D^'^{E,k) = jdE' A^{E\k) ^ ^ ^ 



E' ±ie 

and the limit of free particles is recovered when 

A^{E,k) 5{E^ - El{k))Q{E) 
with the bare "on-shell" energy 

E^{k) = Jk' + ml. 



E + E'±ie 



(13) 
(14) 
(15) 



For the interacting system, the propagator is best expressed through retarded 
and advanced components of the polarization since their analytical structure 
is well defined. One has 



U^'^{E,k) 



Re(n(E, k))TiT^(^{E, k) 
E-E' ±ie 



(16) 



6 



The Schwinger- Dyson equation can be solved directly using this decompo- 
sition. In our formalism it is then diagonal, corresponding to two complex 
conjugate equations for the retarded and advanced propagator. The solution 
is 

ME, k) ^ ^iE,k)eiE) 

^ ' ^ (E2-E2(fc) -Re(n(E,fc)))' + 7rV2(E,fc) ^ ^ 



In the limit of vanishing self energy the free-field result (p^ ) is recovered. 

We now use the fermion propagators @ to obtain the polarization func- 
tion for pions coupled to the fermions. Although we make a one-loop calcula- 
tion, this does by no means imply a perturbative character of the calculation, 
since we can, in principle, use the full Green's function from the fermionic 
sector. All the interaction effects of the baryons are then plugged in their 
spectral functions, missing is only a possible vertex correction term. Such a 
vertex correction term can be parameterized as a medium dependence of the 
coupling constant. 

In this generalized one-loop approximation, one has 



g^''^S^^'\p)g^'^St\p + k)] (18) 



where s = 16/9 is a spin degeneracy factor, k = {E, k) in the integrand and 

gW = = . (19) 

The numerical values for the coupling parameters and masses can be found 
in table 1. 

From the definition (^) one can easily derive the relations and (]TI1|). 
Inserting the fermion propagators given above then leads to a retarded pionic 



JNA 


9' 




Ma 


r 


2 


0.5 


0.14 GeV 0.938 GeV 


1.232 GeV 


0.12 GeV 



Table 1: Coupling constants and masses used in the calculations of this work. 
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polarization function 



6{po - E- zi) + 5{j)Q - Z2) — f20) 

V po- Z2-ie po- E - zi+it) 

+ the same expression with E —>■ —E, e —>■ —e . 
The imaginary part of this expression determines the function a{E, k) as 

'AAiz + E,p + k){nN{z)-nA{z + E)) (21) 
- AAiz - E,p + k) {nN{z) ~ HAiz - E) 



This function rather than (pOD is our starting point for the calculation of 
the full pion propagator, since it is numerically much easier to calculate. To 
include vertex corrections, one could then multiply (t{E, k) by a form factor 
like in ref. |T^, which however we will avoid here. The next step consists in 



calculating the dispersion integral (|TBp. 

Finally, we add another piece which has not been addressed so far. The 
perturbative expression for the pionic polarization tensor in the medium has 
to be corrected for the strong repulsive interaction at short distances. We 
use the phenomenological description 

Note, that in HCNM g' also depends on the system parameters. 



3 Quasistatic zero-width approximation 

A first approximation to the A-hole model is obtained with the free spec- 
tral functions for nucleons and A's, ignoring the finite lifetime of the latter. 



8 



However, since this model is widely used, its critical discussion is necessary 
in this work. Hence we set for the nucleons 



An{E,p) = 5{E-En{p)), (23) 

and neglect the width of the A-particle by using 

A^{E,p) = 5{E-E^{p)) . (24) 

The baryon density of the system is then given by the expression @). 

The quasistatic approximation follows from this by expanding the 6- 
function in the integrand in powers of p/M^. To lowest order, it corresponds 
to the neglection of recoil effects, i.e. by setting 

5{E^{p + k)- En{p) -E)^ 5{E^{k) -Mn-E) (25) 

in the integrand. The 5-functions therefore do not affect the momentum 
integration, and one obtains 



a,,{E,k) = fc'f |/(^ (^^(^^b))-^A(^A(p)) 



X 



[d{E - E^ik) + Mn) - 6iE + E^ik) - M^) j . (26) 

From eqn. (^ follows, that the factor in the curly brackets is equal to 
p% — l/4p^. This finding is in contrast to the equations used in |T0|, |Tl], |1^], 
where the baryon density has been used instead. The difference is crucial for 
the model at higher temperature, where 50-60% of the baryon density are 
A's, hence the factor in the curly brackets is reduced by a factor ^ 2 with 
respect to the baryon density. 

The energy going into the A-production is 



uJA{k) = EA{k) - Mn = ^Jk' + Ml- Mn . (27) 

Note, that in many applications the above expression is expand to lowest 
order in k, thus worsening the asymptotic behavior as ^ oo [0, |Tl|, p^ . 



With this spectral decomposition of the polarization function, the dis- 
persion integral (|I6|) can be performed analytically, giving us the retarded 
polarization tensor in quasistatic zero-width approximation as 




Uf,iE,k) = k^"-[i^] (p^-Va) , (28) 



E^ -ujl{k)+ie' 
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Here e' = sign(i?)e to obtain the proper (retarded) boundary conditions in 
time. The imaginary part is a 5-function, which gives zero contribution to 
the pion propagator. 

The spectral function of the pions therefore consists of isolated poles uj± 
obtained as the solution of 

2 p2.,^ fe'C^A(fc) _ X 



with 



2 



8 ( flVA \ / 1 ^0 



ENAik) = ^ujAik){iUAik)+g'C) . (30) 

Form factors depending only on the momentum k can be absorbed into C. 

For simplicity the momentum arguments are suppressed from now on, 
and thus our solutions are 

1 



^l = ^{El,^ + El± ^{E%^ - ElY + AeCu^ ) . (31) 



The residues of the isolated poles in the propagator are 




En A - E^ 



'{El^-Elf + A^Cu^, 
and the retarded propagator itself reads 



(32) 



DUE,k) 



E - E^^ 



(^2 -ujI + ie') (^2 - cj2 + ie') 

+ .2 \ , • (33) 



^2 - cj2 + ie' ^2 - cj2 + ie' 



For the interpretation of these results we consider the free case, where for 
small momenta E.„ < En a, while at high momenta E^, > En a (see dotted 
lines in fig. 0). 

This can be carried over to the interacting case, where for small momenta 
the spectral weight Z_ > Z+, hence in this case U- is the energy of the pion. 
The other branch in this case is a A-hole excitation. 
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if 


lfc 


< k 




if 




> k 


UJ+ 


if 


\k 


< k 




if 


\k 


> k 



The roles are exchanged, when the momentum is larger. Then, ui^ is 
the energy of the interacting pion with the larger strength Z+, while a;_ is 
associated to the weaker A-hole branch. 

The " crossover" point is determined by the momentum k± where the two 
residues are equal, i.e. 

Elik^) = El^{k^) 

This value is marked in fig. |l] by a vertical line. Hence we can safely set 
pion energy cj^ 

A-hole energy uj^a = { If !^! ^ i!''' • (35) 

The energies uj± are plotted in figure ^ (full lines), for the example value of 
1.69 nuclear density as function of The values of the coupling constants 
and masses used in the calculation are given in table 1. 

Clearly the notion of a dispersion relation with a discontinuity, or jump, 
as implied by the above description, is somewhat unusual and requires further 
comment. We postpone this until we have discussed a way to include the 
width of the branches into the considerations. 

In contrast to our derivation, it has been argued in ref. [|1^, that one 
may introduce the energies of the quasi-pion and the A-hole excitation as 
the two combinations 

^NA = (1 — Z^)uJ- + (1 — Z+)u!-^- = Z+U- + . (36) 

These energies are also plotted in figure 1 (dash-dotted lines) - but their 
choice is unjustified. They do not correspond to any physical excitation of 
the system under consideration. 

To prove this we observe, that from the viewpoint of quantum field theory 
the renormalization of the wave function is done with the factors \/Z±, i.e., 
they are the coefficients of the fields in the dynamical map ||13|. To leading 
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order and for fixed momentum, the dynamical map for the two interacting 
fields can be written explicitly as 

^_ = {zIBb [^^ a] e^v a -4^-Bb \E.\ + . . . . (37) 

Here, ^.^ is the free pion Heisenberg field (thermal doublet, cf. 0), trans- 
formed with the Bogoliubov matrix taken at the proper on-shell energy, ^tva 
is the operator creating a A-hole pair. 

Hence, a Z±-weighted sum of the eigenmode energies ijJ± does not give 
the energy of a physical state, but up to factors \fZ inverts the above trans- 
formation. This can be seen immediately by calculating, with the same jus- 
tification as can be given for eqn. (p6D, the averages of the squared energies. 
The results then are 

Z^Lul + Z+Lul = El 

Z^Lul + Z^Lul = Ej,^, (38) 

i.e., the free energies (dotted lines in figure 1). For all momenta the first 
choice (^) differs only very slightly from the free energies. 

In summary of this comment one can state, that the Z± weighted energy 
average removes the interaction from the system. We therefore feel, that such 
an averaging process should be avoided: the resulting very small modification 
of the pion dispersion relation with respect to the free one can be obtained 
much easier, than by arguments about polarization functions. 



4 Lorentz spectral function 

The use of a free spectral function for the nucleons is justified by the great 
success of the quasiparticle concept in nuclear physics. It can easily be mod- 
ified by introducing a density (and temperature) dependent effective nucleon 
mass, or collective effects due to nucleon-hole excitations [|14|. For simplic- 
ity however we have kept the bare nucleon, and therefore use the spectral 
function as specified in (P^l). 

The A33 resonance in contrast has a vacuum width of 120 MeV. Since 
this is comparable to the pion mass, it should not be neglected in any serious 
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computation. The width of the A spectral function peak also grows with 



momentum |15|, but comparison with scattering data restricts it to about 



twice the above value in the momentum range of interest here |TB 



Even less is known about its medium dependence |]TT|, hence without a 
self-consistent calculation we feel safest to estimate the influence of a constant 
A-width on the pion spectrum. Such an approximation then clearly does not 
have the proper threshold behavior, but can easily be generalized to such 
case. For the following, we therefore use the spectral function 

With the above spectral functions, the actual baryon density instead of (^) 
becomes 

pb = P% + PA 

= 4 y'-^n^(E^(p)) + 16 jdEAA{E,p)nA{E) . (40) 

The imaginary part of the polarization function is 

a{E^k) = J-0^s^A{EM{p) + E,p + k) X 



i^nN{EN{p))-nA{EN{p)+E)j (41) 
— the same expression with E — > —E , 

The complete momentum integral in the second part of (^) can be done 
analytically. For the imaginary part of the polarization tensor in eqn. (p]) 
however, this only holds for the angular integration. Using p = \p\ and 
k = \k\ in the following, we obtain 

a{E,k) = ^J^I^lM^y jp2dp[(nNiE^{p))-nAiE^{p) + E) )IiE,k,p) 

- {nN{EN{p))-nA{EN{p)-E))I{-E,k,p)] (42) 

with the function 

IiE,k,p) = JdnpAAiE + ENip),p + k) 
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-2(En(p) + E) arctan — 

\2x 



(43) 



X = w+ 



X = 



The boundaries to be inserted for x are 



w± = E^{p±k)-EN{p)-E 



(44) 



Remaining for the calculation of cr is therefore the p-integration, which we 
perform numerically. It is followed by the calculation of the dispersion inte- 
gral and by the subsequent correction according to eqn. (|^). The resulting 
spectral function is plotted as function of temperature and energy for five 
different momenta in fig. 2-6. 

Clearly, the spectral function has only one peak (as function of energy) 
at low momenta (fig. 2 and 3) - corresponding to an almost free pion, with a 
momentum dependent width. At momenta in the vicinity of k_[_ according to 
(|3^, this mode is strongly mixed with the A-hole excitation, and two peaks 
appear in the spectral function (cf. fig. 4). 

At higher momenta, this mixing leads to the crossover of the two branches, 
and the diminishing of the peak at lower energies (fig. 5 and 6). It is worth- 
while to note, that the A-hole peak vanishes with increasing temperature at 
all momenta. Hence at temperatures above 0.1 GeV, one cannot speak of the 
A-hole like excitation as a degree of freedom of the system: it is dissolved in 
the medium. 

Our result is, that the location of the peaks corresponds quite well to the 
eigenmodes oj± of the quasistatic zero-width approximation, cf. the dots in 
fig. m. Hence the latter is useful even at finite temperature. 

The same calculations have been performed with a A- width of 2r, which 
therefore allows to study the modifications introduced by a realistic param- 
eterization of the A spectral function in terms of scattering data |]16|. The 



results then are quantitatively, but not qualitatively different from those pre- 
sented in figs. 2-6. They only give an even more pronounced smearing of 
the two spectral peaks at momenta around according to (|3^) (cf. fig. 4). 
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5 Expansion in lowest order F 



The calculation of dispersion integrals with sufficient precision appears to be 
impractical for many purpose. Hence in the following we study a modification 
of the quasistatic zero-width approximation, which takes into account the 
finite width effects in a controlled fashion. 

To this end, we start from eqn. ( p!] ) and perform a quasistatic approxima- 
tion as before by neglecting the p-dependence only in the spectral function 
of the A's. Care has to be taken, however, because of the rather awkward 
factor from the second distribution function. Hence we assume, that Aa is 
still quite peaked around the bare energy value. Then we have approximately 

2^6 / /^a\ / 1 



a,{E, k) = k'- i^p% - -plj {Aa{M^ + E,k) - Aa{Mm - E, k)) . 

(45) 

Note, that still p^^ appears in the coupling factor, rather than as defined 

in (ig). 

Now we follow the approach from and perform an asymptotic expan- 
sion of the A spectral function 

A.(E,p) = S(E-E.ip))-^l.^_^^. (46) 

Here, V denotes the principal value in case one integrates over the energy 
pole. Note, that this equation also holds in case F acquires an E, k, T- 
dependence. Using eqn. (E3), this yields 



TT (E2-^|(fc)) J 

Only the 5-functions contribute to the dispersion integral (0), hence to first 
order in F the polarization function is 
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(before the (yf'-correction according to (p2D). The pion propagator obtained 
with this approximation has two poles in the proper complex energy half 
plane, and is an analytical function in the other half plane. 

Another approximation used frequently is obtained by simply shifting cja 
by iV/2 as in refs. p|, IH: 



However, in contrast to our derivation above, the pion propagator obtained 
with (^) has not the proper structure in the complex energy plane (see the 
comment below and fig. 0). 

After the correction (^21) is performed, the correct pionic spectral function 
obtained with the expression (^91) is, to first order in F, 



AiEk) = ^ ^^"^Q^^) (51) 

Here, to preserve unitarity, the energies in the denominator have been ob- 
tained in slight modification of (0) as 

1 



^'i = 2 (^^A + (r/2)^ + El± ^{E%^ + (r/2)2 - E^y + AeCu;^ 

(52) 

In fig. 1^, we compare the results for the three forms of the asymptotic 
expansion, obtained with eqns. (|48|) (dash-dotted line), (|^) (full line) and 
( |50D (dotted line) to the full calculation from the previous section (dashed 
line). At zero temperature near the point of maximal mixing \k\ ^ k±_ the 
first form ( ^Sf ) is rather crude and hence will be avoided henceforth. The 
approximation (^) however corresponds quite well to the full calculation - 
while the third form (^) contributes unphysical strength at low energies, at 
least for the energy independent F used here. 

While an explicitly energy dependent F might cure this (eqn. (^6]) still 
holds in this case), for our discussion of an energy-independent F we thus 
consider eqn. (|49|) the best choice for an approximate treatment of the A- 
width. For zero temperature, this agreement prevails at all momenta up to 
twice nuclear density, see the upper panel in fig. ^. 

However, the picture changes substantially when looking at finite tem- 
peratures - and this was the main goal of the present work. We find, that 



16 



at finite temperature the asymptotic expansion severely underestimates the 
smearing of the spectrum. This is shown in the small inserts to fig. ^ and 
the lower panel of fig. |^: the discrepancy increases with temperature, but is 
almost independent of density in the range considered here. 

The deviation is connected to the smearing of the Fermi surface, i.e. 
to the quasistatic approximation leading from (|4TD to p5|) rather than to 
the asymptotic equation (|46|) . It can be corrected by fitting an "effective" 
r'(fc, T) to the full calculation. 

In the main body of fig. ^ and the lower panel of fig. ^ we show the 
results of such a fit. Both figures contain a comparison of the full calculation 
from the previous section (thick lines) to the asymptotic expansion (thin 
lines). The asymptotic expansion has been taken from eqn. (PD[), but with a 
modified T'{k, T) - whereas the small inserts to these figures are calculated 
with the r from table 1. 

The effective T'{k,T) is plotted in fig. Of course, such a momentum 
and temperature dependence also arises from a medium-dependent A spectral 
function or its momentum dependent parameterization |T6[. However, we 
find that even when starting from eqn. (|39|) with a constant F, the pion 
propagation feels an effectively medium dependent F'. The changes here are 
of similar magnitude as one would expect from a self-consistent calculation 



of the A-width |]Tl|. The additional width is additive, i.e., if the bare width 
of the A is doubled, the effective width according to our fit does not increase 
by a factor of two. 

It is now clear, how to construct a simple approximation to the pionic 
spectral function in HCNM. This can be used in the calculation of observable 
quantities. However, for the practical purpose of simulations as well as for a 
qualitative understanding, a particle-like picture would be more appeahng. 
We therefore attempt to define an interacting pion propagating in HCNM. 

Such description is far from trivial, since one cannot simply look for poles 
of the pion propagator in the complex plane: they are on an unphysical 
Riemann sheet. Instead, one has to define the particle-like excitations (which 
are not quasi-particles) on the real energy axis. To this end, we follow a 



procedure recently outlined for non-relativistic systems [20|, and write the 



inverse retarded propagator for the pions with real energy as 

kl - El{k) - n«(fco, k) = (A;o - {Sk - ilk)) + i^k + m)) • (53) 
To lowest order, the on-shell energy e and its imaginary part 7 are then given 
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as the solution of the equation 



- ilkf - El{k) - n^(£fc, fc) = . (54) 



Using the asymptotic expansion (^) after the correction according to 
the solution of this equation corresponds to solving two coupled (non-linear) 
polynomial equations. Although this is easily done numerically, one cannot 
guarantee the existence of a unique solution. Indeed, for a broad range of 
momenta two solutions exist. They are shown as full thin lines in fig. ^ the 
shaded areas correspond to the regions Sk + ik > E > Ek — 7fc. 

The fact, that the real part of the energies Sk corresponds quite well to 
the peaks of the full calculation from section 4 proves the validity of the 
approximation (0). Moreover, as can be seen from fig. |Tl], the imaginary 
part of the on-shell energy is comparable to the half-width of the spectral 
function peaks. Thus we conclude, that the definition is an easy way to 
determine the properties of pionic modes in HCNM. 

It remains to clarify, which of the two branches to chose for a particle-like 
pion propagation. A natural choice is, to take as "the pion" that solution 
of eqn. (|5^), which has the smallest imaginary part of the energy. This 
choice is unique, since one of the solutions has an imaginary part increasing 
with momentum, while the other decreases with |fc|. Thus, there exists a 
certain momentum, where one has to switch from one branch to the other 
with increasing or decreasing momentum. As one can find analytically from 
the above equations, this crossover point is given by the momentum value, 
where the two residues in the quasistatic zero-width approximation are equal, 
i.e., by the k±_ defined in eqn. ([3^ ) (cf. also fig. and the vertical lines in 
fig. 0). 

The real and imaginary part of the pion dispersion relation obtained with 



this description are shown in fig. |TT], for several temperatures at 1.69 nuclear 
density. To keep the picture as clean as possible, we have performed these 
calculations with the constant F from table 1. In fig. ^ we show the same 
quantities at three different densities. 

In short, the peculiarity of our definition is a jump in the dispersion rela- 
tion of the interacting pion (To guide the eye, the curves in the figures have 
been drawn continuously.). Such a discontinuity in the dispersion relation 
can be seen as a potential well in momentum space: If the pion momen- 
tum gets smaller than the critical value of k^, the energy drops by a certain 
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amount. This drop diminishes with increasing temperature (see fig. |TT]) and 
decreasing density (see fig. P). The energy released in this process is dis- 
tributed into the medium by the couphng to the broad A-hole excitation (cf. 
fig. 4) . 

The imaginary part of the quasi-pion energy has its maximum at the 
boundary of the well, i.e., at momentum k±. As can be derived, this max- 
imum is max{'-fk) = r/2. The quasi-pion therefore moves through the nu- 
clear medium with a complex velocity, similar to a photon crossing a semi- 
transparent medium. 



6 Conclusions 

In this work, the formalism of Thermo Field Dynamics was applied to the 
A-hole model at finite temperature. The model was analyzed starting from 
a seemingly simple quasistatic zero A-width approximation, which can be 
treated analytically. However, as an important difference to other treatments 
of the model, a diminishing of the coupling was found at higher temperature 
(cf. eqn. (pB]) and its discussion, and compare to refs. |]10|, |12|). 



Furthermore we could show, that the residue-weighted energy average of 
the eigenmodes in this simple approximation does not give a proper descrip- 
tion of pion propagation: one could as well use the free pion in the system. 

We then went to a more elaborate treatment, taking into account a con- 
stant width of the A33 resonance. The numerical calculation of the pionic 
spectral function at finite temperature was performed to establish a refer- 
ence point for simplifications. Of these the most promising was found by 
performing an asymptotic expansion of the A spectral function, leading to 
a polarization function ( ^Q]) that is a straightforward generalization of the 
zero-width approximation. This expansion can also be applied to more so- 
phisticated parametrizations of the A spectral function, where the width is 
E, k, T - dependent. 

Our approximation agrees quite well with the full reference calculation, 
thus supplying a simple expression for the full pion propagator in hot, com- 
pressed nuclear matter (HCNM). Another approximation however was found 
to produce artificial spectral strength at low energies (cf. eqn. (|50|), fig. |^ 
and ref. ||19|| ). 



Since the Bose equihbrium distribution function is strongly peaked at 
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small energies, this strength leads to a large increase of the effective (near-) 
equilibrium occupation number. While such spectral contribution at small 
energies might therefore explain the low momentum pion enhancement in 



relativistic heavy-ion collisions [12, 21 1, it is clearly not a consequence of a 



rigorous treatment of the A-hole model. 

The successful description of phenomena like lepton pair creation in heavy- 
ion collisions requires to use the full pion propagator according to (12), with 
one or the other spectral function as discussed in this work |18|, [19[ . However, 
for cascade-like simulation calculations as well as for other practical purpose, 
a description of "the pion" in terms of spectral functions and matrix valued 
propagators might be impractical. 

We therefore gave a simplified description, i.e. defined particle-like excita- 
tions with a definite (complex) energy. While the real part of this dispersion 
relation coincides very well with the simple quasistatic zero-width approxi- 
mation, our model also gives the spectral width of the pion-like excitation. 
The real part of the dispersion relation obtained in this way includes a discon- 
tinuity, i.e., a potential well in momentum space - and the absorptive part 
of the dispersion relation is maximal at the boundary of the well. Hence, 
this simplified description of pionic modes in HCNM corresponds to a kind 
of "pion-optical" picture. 

One of the implications of this picture is, that the (near-)equilibrium 
occupation number of pions is enhanced for pion momenta just inside the 
potential well, while it is decreased for momenta just above the step (cf. 
eqn. (33)). Further implications, e.g. on pion spectra obtained from nuclear 
collisions, will have to be investigated. 

Furthermore one has to add, that a fully self-consistent treatment of the 
medium- and temperature dependent A-width is desirable. Without such 
considerations, a simple parameterization of our results as function of tem- 
perature and density is meaningless. However, since we derived a quite large 
effect of temperature on the pion width, small changes from such calculation 
will not modify our conclusions. In essence, one has to solve the Fock prob- 
lem for pions, nucleons and A's - and for this, solutions exist that can be 



supplemented by the pionic spectral functions we derived [14 



We did not discuss negative energy states and renormalization problems 
in the present paper. It is nevertheless clear that they play an important 
role for model consistency: Although the A-N continuum begins at much 
higher energies than discussed here, the polarization function has to fulfil 
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the dispersion relation (|T6|). Hence, even at low energies mesonic spectral 
functions are influenced by the continuum |22|. 
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Figure 1: Pionic dispersion relation at pb = 1.69 nuclear density. 
Full thick lines: Quasistatic zero- width approximation, eqn. (2£) 
Dotted lines: Free energies E.,^ and E^/\, vertical lines k± from eqn. p^ ) 
Dash-dotted lines: Z± weighted energies, eqn. (3£) 

Dots: Location of the peaks in the spectral function of the full calculation, (42). 
Full thin lines and hatched area: Asyng^^otic expansion with on-shell definition, 
eqn. (54). Hatched is the region within =b 7^. 



Figure 2: Pionic spectral function at |fe| = 50 MeV, = 1.69 nuclear density. 
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Figure 3: Pionic spectral function at |fc| = 200 MeV. 
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Figure 4: Pionic spectral function at |fc| = 350 MeV. 
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Figure 5: Pionic spectral function at |fc| = 500 MeV. 
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Figure 6: Pionic spectral function at |fc| = 650 MeV. 
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Figure 7: Pionic spectral function at T=0 and pb 
Full line: Asymptotic expansion, eqn. (|^) 
Dashed line: Full calculation from section 4 



1.69 nuclear density 



Dotted line: eqn. (50), Dash-dotted line: eqn. 
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Figure 8: Pionic spectral function at T=0.1 GeV and pb= 1.69 nuclear density 
Thick lines: Full calculation from section 4 
Thin lines: Asymptotic expansion, eqn. (p9|) 

Temperatures 50 MeV (full hues), 100 MeV (dashed) and 150 MeV (dash-dotted) 
Main panel calculated with r'(fc, T) from fig. p!o| . 
Insert calculated with r=0.12 GeV. 30 



Figure 9: Pionic spectral function at T=0 and T=0.1 GeV 
Thick lines: Full calculation from section 4 
Thin lines: Asymptotic expansion, eqn. (pg|) 

/9fe= 1.69 (full lines), 0.92 (dashed) and 0.47 (dash-dotted) nuclear density. 
Lower main panel calculated with r'(fc,r) from fig. [To| . 
Top panel and insert calculated with r=3(D.12 GeV. 



Figure 10: Effective tliermal A-widtli T'{k, T) in GeV at = 1.69 nuclear 
density. 

Both panels show the same data, contour spacing in the lower panel is 10 MeV. 
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Figure 11: Pion dispersion relation at pf, = 1.69 nuclear density. 
Upper panel real part e^, lower panel imaginary part jk of the effective pion energy 
from (g). 

Temperatures 1 MeV (full lines), 50 MeV (dashed), 100 MeV (dash-dotted), 150 

MeV (dash-double-dotted) and 200 MeV (dotted). 

Full thin line: free pion, vertical lines k^^rom eqn. (|3^. 



Figure 12: Pion dispersion relation at T=0 and T=0.1 GeV. 
Upper panel real part Sk, lower panel imaginary part 7^ of the effective pion energy 
from (g). 

Temperatures 1 MeV (thin lines) and 100 MeV (thick lines). 

Baryon densities pb= 1.69 (full lines), 0.92 (dashed) and 0.47 (dash-dotted) nuclear 
density. 34 
Dotted line: free pion 



